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1 Introduction 

This paper presents the major resuh of my doctoral dissertation written at Columbia 
University [11], with Walter Neumann as my thesis adviser. Known uses of the method 
developed, which allows one to conclusively prove that a 3-manifold has a hyperbolic 
structure, include some of David Gabai's [6, 7] recent work and a paper by Chris 
Leininger [10]. Since the determination that M is complete hyperbolic is dependent 
on there being a solution to a set of equations, we shall first review the development 
of these equations. Every orientable complete hyperbolic manifold of finite volume is 
obtained from an ideally triangulated one by Dehn surgery on some of its cusps. This 
fact is documented by Neumann [15], based on a Thurston preprint [18], so we first 
examine A^, a non-compact 3-manifold that is the interior of a compact one whose 
boundary consists of k tori. can be realized as a gluing of n tetrahedra, cxi, . . . , (T„, 
having k vertices after gluing, with a conic neighborhood of each vertex removed [3]. A 
conic neighborhood of the vertex, v, is described as follows. Let f be a vertex and aj a 
tetrahedron that v belongs to. Take the second barycentric subdivision of the edges of 
aj containing v and let Wi, W2 and be the closest vertices to v for these edges with 
respect to this subdivision. See Figure 1. 

Definition 1.1 • L(jj(v) = triangle having vertices wi, W2 and ws as 
above with respect to v and aj 

• L(v) = [Jt, vortex of CTj La-jiy^ 

l<j<n 

L{v) is called the link of v 

• Ucrj(v) = tetrahedron having vertices v, wi, W2 and W3 

• conic neighborhood of v = |J„ vortex ot Ua^iv). 

l<j<n 
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Every vertex is identified with a cusp of and its link is a torus. These truncated 
tetrahedra resulting from the removal of the conic sections can now be treated as ideal 

hyperbolic tetrahedra, so there exists a hyperbolic structure on A^\l-skeleton of N. In 
order for N to have a hyperbolic structure, there must be consistency across the 1- 
skeleton. The conditions for this to happen are embodied in the consistency equations 
and will be described in detail in Section 2, "Identifying the Equations". 

Completeness applies to the cusps. Once a hyperbolic structure is identified, it 
induces a similarity structure (i.e., a (C, AfT(C)) structure) on each of the k tori, 
Ti, . . . ,Tfc. If the similarity structure of a torus identified with a cusp is Euclidean, 
N will be complete at that cusp [3] . This occurs when the image of the holonomy of 
the similarity structure for the torus consists entirely of translations, or equivalently, 
has at least one non-trivial translation [3]. A holonomy of a similarity structure for a 
torus, T, is a map 9 such that 9 : 7ri(T) Aff(C) [3]. The conditions for the image of 
9 to consist entirely of translations are presented by the completeness equations which 
will also be discussed in Section 2, "Identifying The Equations". 

Once we establish the conditions for cusps of N to be complete, we turn our at- 
tention to the manifold M, obtained from N by Dehn surgery on some of the cusps. 
Assume h cusps remain unsurgered, so there are k — h surgered cusps. M must satisfy 
the consistency equations; however, there are now only h cusps that must be shown to 
be complete, so we only need the completeness equations referring to these h cusps. 
The remaining k — h surgered cusps must result from Dehn surgery with co-prime 
coefficients (pj, qi) ior 1 < i < k — h where (pj, qj) and the holonomy of the similarity 
structure of Tj are joined in one equation [3]. 

Once the equations needed to prove a manifold complete hyperbolic are identified, 
we set up the machinery to test whether a solution exists in Section 3, "How to Test for 
a Solution." The method described there concludes the proof of the following theorem, 
which is our main result. 

Theorem 1.2 Let M he a manifold and assume there are n tetrahedra in the triangula- 
tion of M according to SnapPea [21]. There are n equations, {fi{z) = | : C" — >• C} 
for 1 < i < n, whose simultaneous solution will guarantee that M is complete hy- 
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perbolic. If SnapPea finds an approximate geometric solution to these equations, let 
a — (oi, . . . , a„) be an approximate geometric solution generated by SNAP [8] on the 
SnapPea manifold file for M. Let bi = fi{a) for 1 < i < n and f : C" with 

f{z) = {fi{z), . . . , fn{z)), so /(a) = b = (61, . . . ,bn). Then there is L > such that 
there is a genuine solution to the equations, making M complete hyperbolic when the 
following inequality is true: 

< I ^ 

' - 2L|/'(a)-ip 

We devote the final section to examples. Every manifold in the cusped census of 
SnapPea has been examined and the results are reported in Section 4, "Examples." 
However, for detailed discussion, three examples are presented. There are simple ones, 
such as the figure 8 knot complement and Dehn surgery on the Whitehead link comple- 
ment. There is also a complicated link complement with 4 cusps and 32 tetrahedra. In 
uncomplicated cases, it is sometimes possible to show that a knot or link complement 
has a complete hyperbolic structure using means other than the SnapPea approxima- 
tion. Thurston has proven that the figure 8 knot complement has a complete hyperbolic 
structure, and shown when a {p,q) Dehn filling has the same property [20]. Neumann 
and Reid have done the same for Dehn fiUings of the Whitehead link [14]. However, 
when it comes to complicated knots and links, until now, it may have been impossible 
to definitively determine whether this structure exists. For several years Leininger had 
withheld publication of his paper devoted two very large links, one of which is the last 
example [10], because he could not prove that their complements have a complete hy- 
perbolic structure. The paper has now been released using the method presented here. 
So far, every manifold that has an approximate solution with respect to a geometric 
triangulation in SnapPea that has been tested by this method has been verified to have 
a complete hyperbolic structure. 



2 Identifying the Equations 

Let aj be an ideal hyperbohc tetrahedron as described in Section 1, "Introduction", 
and pick an edge e such that Wi & e and prior to truncation, e ended in the vertex v, 
as in Figure 2. Then L„.{v), the triangle with vertices Wi. W2 and naturally has a 
similarity structure as the triangle in C with vertices 0, 1 and z (see Figure 3) [19, 15, 3], 
and the dihedral angle at e will be arg(2;). Clearly, z must be in C+, the upper half 
plane in C. The modulus of L^^ (v) with respect to wi is z, so that the inner angle of 
the triangle at wi is arg(2;). The modulus of aj at edge e is z. The only other modulii 
at the other edges of aj will be either 1 — {1/ z) or 1/(1 — z), so z uniquely describes 
aj in the upper half plane. There are six edges with opposite edges having the same 
modulus [15, 3, 17]. See Figure 4. 
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2.1 Consistency Equations 

In order for N to be hyperbolic, if e is an edge of N, the tetrahedra gluing together at 
e must close up around e. That is, the product of all the edge modulii associated with 
e (different modulus for each tetrahedron e belongs to) must be e^'^*, assuring that the 
sum of the arguments is precisely 2tt. Any of the three distinct edge modulii of a tctra- 

r' // 

hedron, aj, can be expressed as ±-2/(1 — ZjY^ with {r'^^r'-) e {(1,0), (—1, 1), (0, — 1)}, 
so the gluing requirement at edge e is 



where r^- = r'j = if aj does not contain e. A tetrahedron can have more than one edge 
glued at e so r'j and r" can take values between —2 and 2. The Euler characteristic of 
N is zero, so it can be shown that N has n edges [15]. Thus, the n edge equations can 
be expressed as 



±1 



,n). 



(1) 



They are referred to as the consistency equations. The existence of a solution is suffi- 
cient to make N hyperbolic. We rewrite them as log equations because they are easier 
to use this way and it reflects the fact that the sum of the arguments of the modulii 
at each edge is exactly 27i [13]. 



Ci G Z 



(« = l,...,n) 



(2) 



Let R, C and R be the following matrices. 
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Proposition 2.1 // rank H = p, then the space of solutions to the consistency equa- 
tions can be defined by exactly p consistency equations. 

Proof. Let rank R = p < n, so, without loss of generality, we can assume the first 
p rows of R are linearly independent. For s > p, there exist A| e C for 1 < i < p such 
that 



i=l 



i=l 



1=1 
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Assume we have a solution z — {zi, . . . , Zn) to the first p consistency equations. Then 



Thus, 



Hence, 



i=i j=i 



p 



E ((E + (E ^^^^o - ^^o) - (E ^i'^^""' = 0- 

j=l i=l i=l 1=1 

This is the same as 

n 

Therefore, the last n — p consistency equations are determined by the first p, so we 
only need the first p equations to determine hyperbolicity. ■ 
In [15, 3] it is proven that for a complete hyperbolic manifold, rank R = n — 
k. However, we need to prove hyperbolicity. Neumann's work in Combinatorics of 
Triangulations and the Chern-Simons Invariant for Hyperbohc 3-Manifolds [13] tells us, 
without a priori knowledge of hyperbolicity, that rank H = n — k, and C is determined 
by R, so rank H = n — k. This will be explained in Section 2.3, "Matrix Rank". Then, 
by the above proposition, we only need n — k consistency equations to determine 
hyperbolicity. 



2.2 Cusp Conditions 

We now look at the k cusps of N. Details of the following discussion can be found in 
[3]. Let Ti be the torus associated with the i*^ cusp. Select 2 simple oriented loops, 
rrii and Zj, on Tj, representing the 2 generators of the fundamental group of T,. Further- 
more, rrii and li can be chosen as simplicial loops with respect to Tj's triangulation. 
Such a loop is composed of segments where each segment is an edge of some triangle 
Laq{v) C L{v) = Ti, as identified earlier when describing the triangulation of N. Let 7 
be any simple simplicial oriented loop on Tj consisting of d segments, si, . . . ,Sd, and d 
vertices, Wi, . . . , Wd, where Wr is the vertex at the end of Sr as well as at the beginning 
of Sj.+i for 1 < r < ri — 1 and w^, is the vertex at the end of and beginning of Si. 
See Figure 5. We lift 7 to C = M^, the universal cover of Tj, starting at the beginning 
of Si and map it to C by way of the developing map [17, 20]. The resulting curve will 
consist of d straight segments, si, . . . ,Sd, joined at the vertices for 1 < r < (i — 1, 
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Figure 5: Simple Simplicial Loop, 7, on Torus Tj 




Figure 6: Developing Map Image of 7 



as in 7, except at Wd, which does not necessarily connect to the beginning of Si. So 
it starts at the beginning of Si and ends at the end of §4. Repeat the development 
map process, starting at the end of Sd and let s'l be the first segment this time, so 
Wd is the vertex between Sd, and s'l. See Figure 6. Call this curve 7. Aff(C) can be 
regarded as C x C* with (a, 6) G C x C* such that it represents a + bx, an affine map 
of C. The dilation component of (a, b) is b. Thus, if an oriented triangle in C has two 
edges ei and 62 where ei ends in the vertex x, and 62 begins at x, and the modulus 
of the triangle with respect to x is y, then the one and only orientation preserving 
similarity of C that takes Ci to 62 has dilation component equal to —y. Remember, 
the modulus of the triangle with respect to x is defined so that Ci is identified with the 
edge from to 1 and 62 with the edge from to ?/ in the triangle with vertices (0, 1, y). 
If Xri, . . . , Xrp^ are the vertices of the pr triangles, L^,^-^ {v), . . . , L„^^^ {v), that touch 7 
at Wr, as in Figure 5, we get Pr corresponding triangles, Lo-^j (v), . . . , L„^^^ (v), touching 
7 at Wr with Xri, . . . , Xrp^ the respective vertices of these triangles at Wr- The ordering 
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is such that Sr is the first edge of L„^^{v), and s^+i is the second edge of L^^^^{v), 
at Wr unless r = d, and then s'l is the second edge of L^^^^{v). See Figure 6. If the 
corresponding triangle modulii at Wr are Uri-, ■ ■ ■ lUrpri then the dilation component of 
the affine map that takes Sr to s^+i is — 11^=1 Orientation is responsible for the 
"— " in the product. Hence, the affine map that takes si to s'l has dilation component 
of nr=i(~l) nrii?/" = (-l)'^nr=infliZ/"- Note that the modulus of L^^^'") 
for 1 < i < pr is the same as the modulus of Lo-^^(f) at Xj-i for 1 < i < p^-, and this 
latter modulus has already been identified as either Zj^ 1/ (1 — Zj) or 1 — 1/ Zj for some 
1 < J < ^- Therefore, the dilation component of the affine map that takes Si to s'-^ is 
of the form 

n 

±i\{4{l-z^'L 

The holonomy of the (C, Aff{C)) structure on Tj is a map 9 : 7ri(Tj) — > Aff(C) such that 
if [7] is the element of 7ri{Ti) represented by the loop 7, then 9 takes [7] to the affine 
map that takes Si to s[. This is a homomorphism that is well defined up to conjugacy 
class within Aff(C). However, any two elements of Aff(C) within a conjugacy class 
have the same dilation component [3], so the map 

i/ji : ni{Ti) C* such that 

is a well defined homomorphism. ^([7]) will be a translation if its dilation component 
is 1, so ^([7]) will be a translation when ■0i([7]) = 1- 

We now look at loops nii and /j. For simplicity of notation, we also refer to the 
corresponding generators of 7ri(Tj) as mj and li so 

n 

Aim,) = ±ll[z';''^{l-z,)< (3) 

Uh) = ±ifl4^{i-z,y^^ (4) 

If the triangulation of Tj causes to be a simplicial loop with d segments and d 
vertices, then its holonomy will be a non-trivial translation when ipilrrii) = 1 and the 
sum of the arguments of the modulii at the d vertices of rrii is dn [3] . Rewriting in log 
form, these requirements are expressed as 

n 

^{m'ij log{zj) + m'-j log(l - Zj)) = Cmini with Cmi e Z. 
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Similarly, one can identify the log equation which sets the condition for the holonomy 
of li to be a non-trivial translation. It can be expressed as 

n 

^{I'ij log(2j) + I'ij log(l - Zj)) = cuni with cu e Z. 

When the holonomy of the affine structure on Tj has at least one non-trivial translation 
in its image, the afiine structure is Euclidean [3]. But a Euclidean structure on Tj means 
that the i^^ cusp is complete [3], so the completeness equations for all of the k cusps 
are 

n 

^(m^ log{zj) + m"j log(l - Zj)) - c^jTri = {i = l,...,k). (5) 

Now consider a hyperbolic manifold, N, with k cusps where h of the cusps are 
complete, so the above completeness equations hold only ior k — h + 1 < i < k. Let 
Ti be a torus associated with one of the k — h non-complctc cusps. If pi and qi are 
co-prime integers, {pi, qi) Dehn filling can be performed on this cusp. In the literature, 
this process is frequently referred to as Dehn surgery, but it is really a filling. In this 
case, Pifrii + qik is the generator of 7ri(7i) that is killed by Dehn filling. In order to 
extend the hyperbolic structure on N to the Dehn filling at this cusp, we need [15, 3] 

n 

n 

q, ( log(^,) + log(l - Zj)) - Cum) = 2m. (6) 

That is 

n 

({Pim'ij + qil'ij) \og{zj) + ipim'-j + qiQ log(l - zj)) = c^ini with Csi G Z. 
Therefore, if the equations 

n 

{(Pi'^ij + ^i^ij) log(^i) + (Pi^ij + log(l - Zj)) = Csim (7) 

(i = 1, . . . , /c — /i) 

are satisfied, M, the manifold derived from by Dehn filling on the k — h cusps, will 
be hyperbolic near these cusps. 
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The last step in identifying the equations is the selection of the appropriate n — k 

consistency equations. Let s'^ = Piin^j + qil'^j and s'-j = Pim'-j + Qil'lj, and define the 
matrices M, L, S and Mh as R is defined on page 5 so that 



M = 



m 
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We will see that rank U = A;. We can select n — k consistency equations so that their 
rows in R are linearly independent, and when concatenated with U, give an n x (2n) 
matrix of rank n. The reasons for this are a consequence of [13], and will be explained 
in Section 2.3, "Matrix Rank". We will assume, without loss of generality, that the 
last n — k out of n consistency equations are the ones we want. 

In summary, we have 
n — k consistency equations, 

n 

^iKj log(%) + r'-j log(l - Zj)) - CiTii = {i = k + l,...,n), 
k — h surgery equations. 



— CsiTli — 



l,...,A;-/i), 
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and h completeness equations, 

n 

^{T^'ij log(^j) + "^ij log(l - ^j)) - (^miT^i = (i = A; - /i + 1, . . . , /c). 

giving a total of n equations that must have a simultaneous solution to make a manifold 
complete hyperbolic. 

2.3 Matrix Rank 

In [13], Neumann has constructed a chain complex, J, and described its homology. 

Using the terminology of Section 1, "Introduction," with respect to the triangulation 
of N and M, let K be the gluing of the n tetrahedra, (7i, . . . , (T„. The modules of the 
chain complex are Cq, Ci and J, where 

1. Co = Z module generated by the k vertices of K. Each vertex will be associated 
with a cusp of and the torus that is the link of the vertex. 

2. Ci — Tj module generated by E-i, . . . , E^, the n edges of K. 

3. With regard to J, for each tetrahedron, cxj, label the edges as e^i, . . . , e^e according 
to the associated parameters as: 

eji = z ej2 = izi e^-s = 1 - i 

ej4 = z = 33^ 6^6 = 1 - - 

Let JcTj = Z module generated by the six edges of Cj with the relations e^v — 
e-j(T+3) = for 1 < r < 3 and eji + ej2 + ejs = 0. Thus, opposite edges of the 
tetrahedron are represented by the same element of J^-. , and e^a can be defined 
in terms of e^i and ej2- This means that eji and 6^2 generate the Z module, ■ 
Let 

J- II ^ 

l<j<n 

The chain complex sequence is 

J: ^ Co A Ci 4 J ^ Ci ^ Co ^ 0. 
We have a, /3, a* and P* defined as follows: 

1. q; : Co ^ Ci, where a takes a vertex to the sum of the edges containing the 
vertex, with an edge counted twice if both ends of the edge are at the vertex. 
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2. (3 : Ci ^ J can be defined by letting 

Ei^ J2 ^or 

l<j<n i<T<6 

Ej^ is identified with ej^ 

We have the sum J2 i<r<6 ejr G Jaj because more than one edge of aj 

is identified with Cjj. 

can be identified with Ei. 

3. To define j3* : J ^ Ci, note that for each Cj, we have the edge set {cji, . . . , ejg}. 
Let p : {cji, . . . , eje} — >• {Ei, . . . , En} be such that picj-,-) = E^ when e^^ is 
identified with the edge E^. Then, let 

I3*{ejr) = P(ej(r+i)) - p(ej(r+2)) + p(ej(r+4)) - p(ej(r+5)) (indices mod 6) 

That is, 13* takes e^v to the alternating sum of the edges of N identified with the 
edges of Uj that touch Sjr- 

4. a* : Ci ^ Co, where a* sends an edge, Ei, to the sum of its end points. 

N is the interior of a compact manifold, A^, whose boundary is the union of the k tori, 
Ti, . . . , Tfc, that are the links of the vertices of K. 

Lemma 2.2 When tensored with Q, the sequence, J , is exact except in the middle, 
where its homology is Hi{dN;Q) = lli<j<fe -f^il^^ii Q)- 

For a proof, see [13]. We use this to compute the rank of R. However, we will use the 
original chain with coefficients in Z to show that the rank of the matrix obtained by 
concatenating U, as defined on page 10, with n — k linearly independent rows of R, is 
n. 

2.3.1 Rank of R 

The matrix of the hnear transformation, /3, is closely related to R*, the transpose of 
R, and they have the same rank. Since rankR = rankR*, rankR = rank of the 

matrix of /3. The edges . . . , basis of Ci vector space, so the vectors 

f3{Ei) for 1 < i < n are the columns of the matrix of jJ. From the definition oi jJ, we 
see that in J^., 

P{Ei) = ^ Cjr modulo relations on J. (8) 

1<T<6 

is identified with ej^ 

Thus, if: 

Cji or occur, it means E^ is identified with the Zj parameter 
ej2 or occur, it means Ej, is identified with the parameter 
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Cjs or CjQ occur, it means Ei is identified with the 1 — ^7 parameter. 



In Jaj, ejs = -eji - 6^2; also, 1 - 



Hence, the sum of the coefficients of 



Cji in (3(Ei) is r'^j, the sum of the exponents of zj with respect to the edge Ei in the 
consistency equations, and the sum of the coefficients of 6^2 in (3{Ei) is —r'/j, which 
is —1 times the sum of the exponents of 1 — Zj with respect to the edge Ei in the 
consistency equations, as seen on page 5. Consequently, R, the 2n x n matrix of /? is 



R 



11 

' 11 
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In 
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We see that rankR = rankR*, so rankR = rankR. By definition, the rank of R 
is equal to the dimension of the image of (3. By Lemma 2.2, a is injective, making 
dimim(Q!) = dim Co = /c, and im(a) = ker(/3). The matrix of /3 this way would still 
be R, so 



rankR = dimim(/?) 

= dim Ci — dimkernel(/9) 

= dimCi — dimim(Q;) 

= 11 — k 



Therefore, rankR — n — k. Let R, C and R be the matrices associated with the 
consistency equations, as on page 5. Consider the matrix equation R-x = — C. In [13] 
it is proved that there is an x G Q^" that is a solution. Then — C is a linear combination 
of the columns of R, so R concatenated with — C has the same rank as R since row 
rank is the same as column rank. That is, 

n — k — rankR = column rankR = column rank (R|— C) 
= column rank (R|C) = column rankR = rankR, 

so rankR — n — k. Let R/3 = matrix consisting oin — k linearly independent rows of 
R. 



2.3.2 Rank of (S|Mh|R/3) 

For now we will include all k cusps of N. Let Si{dN) = Z module of simplicial 1- 

chains, Zi{dN) = Z module of 1-cycles and Bi{dN) = Z module of 1-boundaries. Let 
Cjr G Jcr^ for r = 1, 2. If the two vertices at the ends of ej^ in (Tj are t'^vi and Vjr2, let 
(jTi and CjT-2 be the respective edges of L^. (f jri) and L^. {vjT-2) that do not intersect 
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Cjr- Do the same for ej{r-\-z)-i so we have four 1-simphces identified in dN . They are 
CjTi) 0r2, 0(r+3)i ^-nd Cj(r+3)2, with One for each vertex of Uj. Now define 70. 

70 ^ ^i(9iV) 

Cjr Cjrl + Cjr2 + Cj(r+3)1 + Ci(T+3)2 

We have, by [13], 

70: im(/3) ^ Bi{dN) 
70 : ker(/3*) ^ Z,{dN) 

so there is the induced map 

7:ker(/?*)/im(/?)^i7i(a]V)= ]J //i(T,) 

i<i<fc 

Next, let 60 : Hi{dN) — > J be defined as follows. Let F be a simple simplicial loop 
on the torus, Tj, associated with the i^^ cusp of A'^. In figure 5, 7 is such a loop. 
Each vertex, Wr, of 7, is the vertex of Pr triangles L^^^{v), . . . , L(r^^^{v) where Tj is 
the link of v, a vertex of K. Define the simple cellular path F, by starting at the 
midpoint of the edge of Lo-ii(f) that ends in Wi but is not Si. Continue across the 
{Lo-j {v)}2<q<p^-i by crossing from one triangle to another at the midpoint of the 
edges that have wi as a vertex, ending at the edge of L„j^^^{v) that is not S2- Then 
continue across L(^^p^(v) = L^^-^{v) to the edge of La-^-^{v) that has W2 as a vertex 
but is not S2- Repeat the process until the loop is closed by going from the edge of 
Lfj^^^{v) = Lo-ii(f) that contains but is not Si to the starting point. When F crosses 
(v) for 2 < q < Pr — 1, it goes counterclockwise around the vertex vertex of 

La^^{v), and when it crosses L„^^^{v) = L(^^^^^^^{v), it goes clockwise around the vertex 
of this triangle that is opposite to Sr+i- When one of these vertices belongs to the 
triangle La-^^{v), the vertex is associated with an edge, Crqr, of cr^g for some 1 < r < 6, 
as defined at the beginning of Section 2.3, "Matrix Rank", and this edge is an element 
of J„^^ C J. To each of these edges assign a "+" if F goes around its corresponding 
vertex counterclockwise, and a "— " if F goes around its corresponding vertex clockwise. 
F is homotopic to F, so we can define So : Zi{dN) — > J such that 5o{T) = 5o{T) is the 
signed sum of these edges in J. That is, 

^o(r) = Yl ^-^y^r,r (9) 

l<r<d 
2<q<j>r 

where t — Q when e^^r is assigned a "+" and t — 1 when Cj-qr is assigned a "— ". In 
J(jrqi ^rqr — ^rq{T+3) for 1 < T < 3 with the last subscript mod 6, and e^gi + €^^2 + erg3 = 
0, so —Crgr — erq(T+i) + Grq{T+2) with the last two subscripts mod 6. Therefore, when 
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Crqr IS assigned a "— " , we substitute erq(T+i) + erq(r+2) with both subscripts mod 6. 
Hence, 



l<r<d 
1 <g<Pr 



where e^qr is an edge of a^g that is associated with Wr-, a vertex of L„^^{v) and 
is a vertex of the simple simphcial loop F in Tj. The relations of J also mean that 

5o(r)= J] g'-ej, + 9';^ej2 (11) 

_ l<J<n 

r crosses Lcrj (v) 

where, with respect to aj, g'.r^ is the number of occurrences of the Zj parameter minus 
the number of occurrences of the 1 — ^ parameter and g"^ is the number of occurrences 
of the Tz— parameter minus the number of occurrences of the 1 — - parameter in 

1 Zj Zj 

Equation 10. 

Now let nii and li for 1 < i < A; be the meridianal and longitudinal simple simplicial 
loops on Tj, as in Section 2.2, "Cusp Conditions". We get corresponding rfii and Ij, 
constructed as F was, where and U are homologous to rrii and Ij, respectively. So 
rrii and li are the generators of Hi{Ti) and their image under 5q are two columns of 
ML, the matrix of Sq. These two columns are of the form 

9mi = {g'lrui , dirui ^ ■ ■ ■ ^ 9nmv 9nmi ) with g'^^. = and g'~. = -m'-j 
9k = (4'<I,'---'^n,,<i,) with 4 = and 4 = -Zf^. 

where m^pm'lj and I'ljil'lj are the components of the matrices M and L from "Cusp 
Conditions" on page 10. Let ML be M concatenated with L. For each generator 
of Hi{dN) — lli<i<fe -f^i(^j)> there is a column in the matrix of Sq, so ML has 2k 
columns and 2n rows, where the (2j — 1)*^ row of ML is equal to the j*'' column 
of ML and the 2j*^ row of ML is (—1) times the {n + j)*^ column of ML. Thus, 
rank ML = rank ML * = rank ML. The next step is to show that rank ML = 2k. We 
have im((5o) C ker(/5*), with So{Bi{dN)) C im(/5), so there is the induced map 

6 : Hi{dN) ^ ker(/5*)/im(/5) 

Now 7(5 : Hi{dN) Hi{dN) is multiplication by 2 [13], so 6q must be injective. 
Consequently, the matrix of has maximal rank, which is 2k, making the 2k vectors, 
{9mi, 9ii}i<i<k, linearly independent. 

M is derived from by the Dehn filling oik — h cusps of with fiUing coefficients 
of {pi, Qi) foT 1 <i <k - h. Let ^5. = pi^*^. + Qigi^ for 1 < i < k - h. 

Lemma 2.3 The k + h vectors 

r— * — *— * — *— * — 

\9si 7 • • • 7 9sk-h ' 9mk-h+i ' • • • ' 9mk 1 9lk-h+i ' • • • ' 

are linearly independent. 
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Proof. Assume otherwise. Then there exists for 1 < i < A; — /i and and ipu 
iork — h + l<i<k such that 

l<i<k-h k-h+l<i<k 

= X] (l^siiPigrrii + qSi) + iCmigmi + vugii) 

l<i<k-h k-h+l<i<k 

= ^siPigmi+ ^ <l>siqigii+ iCmigmi + vugh) 

l<i<k-h l<i<k-h k-h+l<i<k 

We have just seen that {gm^, gii}i<i<k, is hnearly independent, so = V«i = 
k — h+l<i<k and = (psiQi = for 1 < i < k — h. But at least one of Pi or 
is not 0, so (psi = ioT 1 < i < k — h. ■ 
Since rankR — n — k, select n — A; linearly independent vectors in im(/3) that are 
columns of the matrix R, and denote them by gpi ior k -\- 1 < i < n. Observe that 
im(5o) n im(/9) = {0}, because otherwise, there is a non-trivial x G Hi{dN) such that 
%6o{x) = 7o(element of im(^)) G Bi{dN). Then ^6{x) = G Hi{dN). But 7^ is 
multiplication by 2 on Hi{dN), so a; = 0, which is a contradiction. 

Lemma 2.4 Let 

1) S= the 2n X {k — h) matrix whose columns are the vectors gs-, for 1 < i < k — h 
^) Mh — 2n X h matrix whose columns are the linearly independent vectors 
for k — h+l<i<k 

3) = the 2n X {n — k) matrix whose columns are the linearly independent vectors 

gp^, fork+l<i<n 

Concatenate these matrices to get the 2n x n matrix F = (S|Mjj|R^). Then RankF — 
n. 

Proof. Assume otherwise. Then the vectors that are the columns of F are not 
linearly independent, so there are isi ^or 1 < i < k — h, ipmi ior k — h + 1 < i < k and 
(ppi ior k + 1 < i < n, where not all are zero, such that 

= Yl ^''3si + Yl frnigrrii + Yl ^f^'^Pi (12) 
\<i<k-h k-h+l<i<k k+l<i<n 

Therefore, 
= 7o(0) 

= 7o( Y ^^i9si+ Y frnigmi+ Y ^P^dp^j 
l<i<k-h k-h+l<i<k k+l<i<n 

= 7o( Y ^^i9si+ Y V'rniffmi) +7o( Y ^Pidp^j 
l<i<k-h k-h+l<i<k k+l<i<n 

= 7opo( Y ^siiPim + qJi) + Y 'yi'mi^^i)] +7o (element in im(/?)) 

l<i<k-h k-h+l<i<k 
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But 7o(im(/3)) C Bi{dN), so lS(^J2i<i<k-h^si{Pim + qik) +T,k-h+i<i<kVrnim^^ = 0. 

Therefore, Y.i<i<k-h^siiPimi + qik) + J2k-h+i<i<k^rnirrii ^ since 7*5 is injective. 
Hence, 

= 5o(0) 

l<i<k-h k-h+l<i<k 
l<i<k-h k-h+l<i<k 
l<i<k-h k-h+l<i<k 

By Lemma 2.3, for 1 < i < k — h and for k — h+l<i<k are all zero. 
Then, Equation 12 becomes = Sfe+i<i<n However, the ^/j. for /c + 1 < i < n 

were selected to be hnearly independent, so (f)f3i — ior k + 1 < i < n. This is a 
contradiction. ■ 

Corollary 2.5 Each column ofR^^ has a corresponding row in R, the matrix associated 
with the consistency equations. Let Hp be the matrix comprised of only these n — k 
rows of a and let 

( ' 

Then rankF — n. 

Proof. As before, every (2j — I)*'* row of F is equal to the j*^ column of F and 
every 2j*'* row of F is (-1) times the (n + j)*^ column of F. Thus, 

rankF = rankF* = rankF = n 

That is, rank ¥ — n. ■ 

3 How to Test for a Solution 

Let 

n 

fi{zi, ...,Zn) = {^i'^'ij + ^i^'ij) ^Og(^j) + (Pi'^ij + ^i^ij) ^Og(l " ^j)) 

—CsiTii {i = 1, . . . , k — h) 

n 

fi(zi, ...,Zn) = ^i'^ij log(^j) + m"j l0g(l - Zj)) - Cmi^i 
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fi{zi,...,Zn) = ^{r'ijlog{zj)+r'^jlog{l-Zj))-CiTTi 



(i = A; + 1, . . . ,n) 



and let 



/:C" ^ C" such that 

Z ^ (zi,...,Zn) f(z) ^ (fi{z),...,fn(z)). 



Then let 



t'ij = Pim'ij + Qil-j t'lj = p^m'lj + qil'lj t'l' = c^i {i = 1, . . . , k - h) 

tij = ''TT-ij tij = J^iy ti = c^i (i = A; — /i + 1, . . . , /c) 

^^i ~ ^^'i ~ t'i' — {i — k + 1, . . . ,n). 

The resulting components of / are 

n 

Mz,, ...,zn) = 5^(4 log(2;,) + t:;. log(l - zjj) - t^'m {t = l,...,n). (13) 
Then = !k _ for 1 < i < n, so 



dZj 



= P-T^,...,^-T^). (14) 



Let if = C!f:, the upper half plane in C". H is open in C". Each /j is holomorphic on 
H, so / is holomorphic on H [16] . Thus / is smooth on H, with the derivative of / at 

f'{z) = [22], being well defined on H. Since we are only working with 

manifolds where SnapPea finds an approximate solution to / in C" , there is an a e C" 
such that /(a) = b and b is extremely close to e C". We know that det/'(a) 7^ [4], 
so rank f'{a) — n and / is regular at a. Then f'{a)~^ exists. Let 

5 : C" C" 

V ^ \f'ia)-'-v\. 

Since 6, as a function of v, is a continuous function on C", it will attain a maximum 
and minimum on the compact set € C" : |v| = 1}. 



3.1 Kantorovich 

The Kantorovich Theorem [9] provides a test for the solution of /. The relevance of this 
theorem to the solution of / was brought to our attention by Joan Birman after another 
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test had been developed by us. We thank her for teUing us about it. The Kantorovich 

Theorem is usable in our situation because we can identify the quantities used, though 
this is not the case for all functions. The test provides a sufficient condition for a 
manifold to have a complete hyperbolic structure. Consequently, it is possible for a 
manifold to not satisfy the condition and still be complete hyperbolic. 

Theorem 3.1 (Kantorovich) Let U he an open neighborhood of a point, a, in C" 
and f ■ U — > C" a holomorphic mapping with invertible derivative f'{a) at a. Let 
hh — —f'{a)'^f{a), a — a + hh and Uq — B\hh\{5,), the open ball of radius \hh\ about 
a. IfUodU and 

1. The derivative f\z) satisfies the Lipschitz Condition on Uq, with Lipschitz Ratio, 
L 

2. |/(a)||/'(a)-T^<i 

then f{z) — has a unique solution in Uq. 

The Kantorovich Theorem applied to our function, /, works as follows. Let U — H. 
Given a, an approximate solution to f{z) — 0, apply Newton's method to / at a to 
get an even better approximate solution, a. That is, let hh = —f'{a)^^ ■ f{a) and a = 
a-\-hh = (oi + hhi, . . . , On + hkn) so dj = Oj + hhj. Then see if a Lipschitz Ratio, denoted 
by L, can be identified for z G i?|/i/j|(a) so that f'{z) satisfies the Lipschitz condition on 
Uq with L. One way to do this is to find an upper bound, Cijk, on the second partials, 

\didjfk{z)\ for 1 < i,j,k < n ior z e S|M|(a), and let L = y^Ei<ij,fc<„(ciifc)^ [9]- 
This works for us, but in general, the major stumbling block to using this theorem is 
the difficulty in finding this L. Here, |/'(a)~^|, the norm of f'{a)~^, can be either the 
supremum norm, which we will denote by |/'(a)~^|sup, or the length norm, referred to 
as I /'(a) lien, where 

|/'(a)"^|sup = sup|^|=i|/'(a)"^ • v\ 
and if a component of f'{a)~^ is denoted by hij, 

l/'(«)"'|len= J E l^^^l'- 
y l<ij<n 

Now substitute values in the inequality found in the second part of the Kantorovich 
Theorem and see if they pass the test. If so, there is a solution in B\hh\{d). 

3.1.1 Calculate \f'{a)-'^\ 

Supremum Norm:|/'(a)~"'^|sup Let 

B = {f'{a)-^ • : = 1} = e C" : \ f'{a) ■ w\ = 1}. 
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We look at the continuous real valued function fi on the compact set B such that 

fx: B ^ R 
w — >• 

Let S = {v & C"- : \v\ = 1} . Then fi attains a maximum at some w ^ B and the 
function 6 will attain a maximum at some v & S where w — f'{a)~^ ■ v. Now let 
A = f'{a). This is a complex matrix, so 

\Aw\'^ — {AwY{Aw) A — conjugate of A and t — transpose of A 
= {w'A')(A^) 
= w\A'A)w. 

Let D = (A^A). This is a self adjoint matrix so it has real eigenvalues [5]. Then, 

B ^ {w: \Aw\ = 1} 
= {w : \Aw\'^ = 1} 
= {w : w^Dw ^ 1}. 

Using the Lagrange multiplier method to maximize n on B [5], let 

H{'Wi, . . . ,Wn, X) — — A(w*Dw — 1) 



n n 



1=1 1=1 j=i 

In order to find a critical point for i7, all partials with respect to wi, . . . , and 
A must be 0. We set 

dH " ^ 

= -K— ^Wj- A(y^ dijW]) = 



so. 



Q — w — \Dw 
= (I-XD)w 

= {jI-D)w. 
Then Dw — ^, making j an eigenvalue of D. Also, 

= w^Dw - 1. 
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Thus, w*Dw — 1, and substituting j for Dw from above, we have w^j — 1. That 
is, w^w = A. But w^w — \w\'^, so maxji/;!^ = maxA such that j- is an eigenvalue 



of D. Then, 

maxonfiiwr = 



1 



smallest eigenvalue of D 
1 

smallest eigenvalue of A^A 
1 



smallest eigenvalue of f'{aYf'{a) 
By definition, |/'(a)~^|sup = niaXonB|w|, so 

l/'(«)-^up = ^ (15) 

a/ smallest eigenvalue of f'{ayf'{a) 



We calculate the eigenvalues of f'{ayf'{o) using its characteristic polynomial and 
take the square root of the smallest one to get |/'(a)~-'^|sup- 

Length Norm:|/'(a)~^|ien Let the components of /'(a)~^ be (/iy)i<i,j<n- Then 



l<ij<n 



3.1.2 Calculate c, 



Let z G -B|/j/i|(a). Then — a| < \hh\^ so — 5^1 < \hh\^ where Zj — aj = Zj — {aj + hhj) 
since dj = aj + hhj. Figure 7 shows the situation for each j. There are three tests that 
need to be performed before we test for the inequality in the Kantorovich Theorem. 
The entire process stops and Kantorovich tells us nothing about a manifold when any 
of these tests fail. 

Test 1 We want a solution in H, so we require that Im(aj) > |hh|. Otherwise, there 
are z e B\hh\{0') that have l'ai{zj) < 0, and the solution could be one of these z. 

Test 2 

\{zj - ttj) - hhj\ = \zj - {ttj + hhj)\ = \zj - aj\ < \hh\. 

Using triangle inequalities, 

\zj - aj\ - \hhj\ < \{zj - Qj) - hhj\. 

Therefore,! —Oj I — I < \hh\, giving \zj — aj\ < \hhj\ + \hh\. But \hhj\ < \hh\, 
so \zj — Qjl < 2\hh\. Now 

\zj\ = \aj + {zj - aj)\ > \aj\ - \{zj - aj)\. 
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Figure 7: Disc of radius \hhj\ about dj 



Thus, 

\zj\ > \aj\ - 2\hh\. (17) 

We need \aj\ — 2\hh\ > in order to define L, so the second test is to check 
that |hh| < ||aj|. Then, 

Jz'l ^ \aj\ - 2\hh\' ^^^-^ 

Test 3 We do a similar process as in the previous test. We aheady know that jzj— a^l < 
2\hh\. Only now, we use 1 — Zj instead of Zj, so 



|1 - Zjl = |1 - % - {zj - %)| > |1 - % 



l^j -^JI- 



Hence, 

\1 - Zj\ > \1 - aj\ - 2\hh\. (19) 
We need |1 — aj\ — 2\hh\ > 0; a third test is to see if |hh| < ||1 — aj|. Then, 



|1 — |1 — aj| — 2\hh\ 



Remainder of Calculation 



Wc are now ready to look at the second partials. By page 18 we see that for 

z e Bihh\{a), 



Therefore, 



dz, z, 1-z, 



dkdjMz) = for A; ^ J (21) 
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Consequently, 



mM^)\<!H 



\t"-\ 
I I 



\l-z,\^ 



Combining this with Equations 18 and 20 yields 

Using this, Qjfc for I < i, j, k < n is defined as 
Ciik = foij^k 



\f" I 









(1% 




-2\hh\y 



The Lipschitz Ratio, L, can now be identified as 

l<ij<n 



It"- 1 



l<ij<n 









(1% 


-2\hh\)^ 



+ 



/" I 



(23) 

(24) 

(25) 



(26) 



The theorem can finally be applied, testing to see if \f{a)\\f'{a) < |. Since 
b — f{a), this can be rewritten as 

\b\< 



2\f'{a)-^'L 



We really have two tests, one using the supremum norm and the other using the length 
norm. This completes the last part of the proof of Theorem 1.2. 



4 Examples 

The methods presented are implemented by the use of two programs: SNAP to get 
information about the manifold and Pari-Gp [2] to do calculations. We use Pari-Gp 
instead of Mathematica because of its high level of precision. Page 24 begins the 
template for a program that is written in an edit file and then copied into Pari-Gp 
for execution. The template needs to be adjusted for information gotten from SNAP. 
Assume we have a manifold file in SNAP for the manifold, M. Once SNAP is open, read 
in the file for processing. The "pr sol" command will print the type of solution SNAP 
has found. A geometric solution means that the solution is in H. Any other response 
is useless here, so there is no need to go any further. Assuming it is geometric, proceed 
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with setting up the template. Issue the "pr sh" command. SNAP will return the 
transpose of a vector representing an approximate solution to our set of n equations 
for M. The number of components of the vector will be equal to n, the number of 
tetrahedra in the triangulation. Then copy this vector from SNAP to the template, 
replacing [oi, . . . , a„], so that a now has the value of our approximate solution. The 
tilde at the end of the Snap response must be eliminated so that a appears as a 1 x n 
matrix. Next comes the "pr fill" command. SNAP will display a (n+k) x (2n+l) matrix 
where the components of each row are the coefficients of a cusp or consistency equation. 
Assuming M is the result of Dehn filling on h out of k cusps, the first k—h rows represent 
the cusp surgery equations, the next h rows are the meridian completeness equations 
for the unsurgered cusps, and the last n rows are all the consistency equations before 
any have been eliminated. If all cusps are unsurgered, h — k, so the first k rows are all 
meridianal completeness equations. This command to print filling equations is closely 
related to the "pr gl" command which prints the gluing equations. This latter display 
presents the k meridianal followed by the k longitudinal completeness equations for the 
original k cusps before any surgery, and then all of the n consistency equations. But 
it is simpler to use the filling equations, even for manifolds where no surgery has been 
done. Copy this matrix from Snap to the template, initializing the matrix FG. The 
script will then create the matrices F and G, where F consists of the first k rows of 
FG and G consists of the last n rows of FG. The rows of F are hnearly independent 
and the program selects n — k rows from G so that when added to F, the resulting 
matrix has rank n. 

The only further adjustments may be Pari-Gp punctuation to reflect line contin- 
uation. In order to tell Pari-Gp to ignore an end of line from the text editor, a "\" 
followed immediately by Return must end that text line. This is needed with a large 
vector or matrix, so it will probably be needed once the values for a, F and G are 
copied into the template. 



4.1 Template 

/* set precision to 60 (or higher for very large manifolds) from the default of */ 

\p 60 

/* read the file FILENAME into SNAP */ 
/* see that there are h unsurgered cusps */ 

/* (,) (,) (,) (,) */ 

/* print shapes - the triangulation has n tetrahedra */ 

/* enter the shapes as a vector, so it is regarded as a 1 x n matrix */ 

a = 

[ai,...,a„] 
/* find n, the number of tetrahedra */ 
n = matsize(a) [2] 

/* print fiUing equations and use this to initiahze the matrix FG. The first \ 
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k — h equations are cusp surgery equations, followed by the h meridianal \ 
completeness equations, and finally all of the n consistency equations. */ 
FG = 

[Xii , . . . , a;i(2n+l) ; • • • ; , • • • , X(^n+k}{2n+l)] 

/* find n + k, the number of equations derived from the Snap command \ 
"pr fill" */ 

numallcq = matsize(FG)[l] 

/* find total number of cusps, k */ 

k = numalleq - n 

/* initialize F, the cusp equations matrix, using the first k equations from \ 
FG */ 

F = matrix(k,2*(n) +l,n,jj,FG[n,jj]) 

/* initialize G, the matrix of all consistency equations, using the last n \ 

equations from FG */ 
G = matrix(n,2*(n) +l,n,jj,FG[k+n,jj]) 

/* define matrix H by eliminating the last column of F representing the \ 

ni coefficient */ 
H = matrix (k,2*n,i,j,F[i,j]) 

/* define matrix K by eliminating the last column of G representing the \ 
ni coefficient */ 
K = matrix (n,2*n,i,j,G[i,j]) 

/* redefine F and H by adding rows to them from G and K respectively \ 

until the rank of F and H are both n */ 
r = 1 

v(r) = vector( (2*n)+l,l,G[r,l] ) 
t(r) = vector( 2*n,l,K[r,l] ) 

while( n — matrank(H) && (n+1— r), if( (matrank(concat(F,v(r))) \ 

— matrank(F)) && (matrank(concat(H,t(r))) — matrank(H)), \ 
(F = concat(F,v(r))) && (H = concat(H,t(r))), r=r+l)) 

eval(F) 
eval(H) 

/* set up the filling equations as log functions evaluated at a */ 
f(i) = sum( j = 1, n, F[i,j]*log(a[j]) ) + ( j = 1, n, F[i,n+j]*log(l-aLi]) ) \ 

+ F[i,(2*n)+l]*Pi*I 

/* define the vector in C" */ 
b = vector( n, i, f(i) ) 
/* identify the norm of 6 */ 
normb = sqrt( norml2(b) ) 

/* identify A, the derivative matrix for / at a */ 
g(i,j) = ( F[i,j]/a[j] ) - ( F[i,n+j]/(l-a[j]) ) 
A = matrix( n, n, i, j, g(i,j) ) 
/* check that determinant of A is not zero */ 
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matdet(A) 



/* KANTOROVICH PROCESSING */ 
/* find eigenvalues for D = (transpose of A) * (conjugate of A) */ 
D = mattranspose(A)*conj(A) 
wapprox = polroots( charpoly(D,x) ) 
w = real( wapprox ) 

/* cliange b into a matrix to do matrix multiplication */ 
B = matrix(n,l,j,i,b|j]) 

/* define the vector hh and find its length, normhh */ 
hhh = -(A)^(-1)*(B) 
hh = vector (n, j, hhh[j,l]) 
normhh = sqrt (norml2 (hh) ) 

/* perform the first three tests to see if this method is applicable */ 
atilde = a + hh 

/* test 1 to see if fat solution; if j > n */ 
for (j = 1, n, if(normhh < imag(atilde[j]), , \ 

error ("failure at atilde[", j, "]"))) 
/* test 2 to see if Cijj can be defined */ 
for(j = 1, n, if(normhh < (l/2)*abs(a[j]), , \ 

error ("failure at atilde[", j, "]"))) 
/* test 3; other test to sec if c^j can be defined */ 
for(j = 1, n, if(normhh < (l/2)*abs(l — a[j]), , \ 

error ("failure at atilde[", j, "]"))) 
/* identify the Lipschitz ratio, Lips */ 
c(i,j) = (abs(F[i,j])/(abs(ati]) - 2*normhh)~2) \ 

+ (abs(F[ij+n])/(abs(l-atj]) - 2*normhh)~2) 
Lips = sqrt( sum( j = 1, n, sum(i = 1, n, c(i,j)"2) ) ) 

/* identify normAinv, the norm of A"(— 1), using the definition of matrix \ 

norm as the supremum of A" {—l)v for v on the n-sphere */ 
normAinv = l/sqrt( vecmin(w) ) 
/* do the Kantorovich tests */ 

/* find the value that the norm of b must be less than or equal to with \ 

respect to the supremum norm */ 
1/(2 * (normAinv) "^2 * Lips) 
normb <= 1/(2 * (normAinv) "2 * Lips) 

/* find the length norm and the value that the norm of b must be less than\ 

or equal to with respect to the length norm */ 
sqrt (norml2 ( A " ( — 1 ) ) ) 
1/(2 * norml2(A-(-l)) * Lips) 

normb <= 1/(2 * norml2(A~(-l)) * Lips) ■ 
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4.2 Using the Template 

The template is ready to be used. If you want a copy of what has happened, first turn 
the log on in Pari-Gp by typing "\1 logfilename" . Then copy the adjusted template to 
Pari-Gp, wait for the run to complete, and open the log file to see the results. Make 
sure that there are no error messages from the qualification tests described. If there are, 
any further results are of no value. If there are no error messages, a response of "1" to 
either of the Kantorovich inequalities will indicate the manifold is complete hyperbolic. 
A copy of the template can be found at [12]. We now look at three examples. Each 
example will have two sets of data. The first comes from SNAP and the second is the 
result of calculations in Pari-Gp. The vectors and matrices are printed as they appear 
in SNAP. When one of them extends beyond one line, it is edited once copied into the 
template to add the hue continuation character, "\," after each hue before its end. The 
Pari-Gp data has been shortened to 40 decimal places from the calculated precision of 
60 decimal places so as to fit on one line since in these examples, it has no effect on 
understanding the results. 

1. FIGURE 8 KNOT COMPLEMENT 

The simplest is the figure 8 knot complement. We know [20] that this is complete 
hyperbolic already. However, only sufficiency conditions have been presented 
here, so it is nice to see that a manifold we know to be complete hyperbolic does 
not fail the test. 

QUANTITIES FROM SNAP 
n = 2 

a = 

[0.5000000000000000000000000000+0.8660254037844386467637231707*1, 

0.5000000000000000000000000000+0.8660254037844386467637231707*1] 

F=[l, 0, 0, 1, 0] 

^=[2, -1,-1, 2, 0; -2, 1, 1, -2, 0] 

Pari-Gp CALCULATIONS 

|6| = 1.296666384352891444530724934775173278518£;- 28 
L = 4.472135954999579392818347339211785668123 

|/'(a)"^|sup = 1.592226038754547070932399593119376104348 
|/'(a)"^|ien = 1.632993161855452065464856049716587347937 
2\f,(„yH2 L = 0.04410070808503045666350407221846082500302 

.ill ' = 0.04192627457812105680767200627679720162466 

W lien-'" 

2. (9872, 11111) DEHN SURGERY: WHITEHEAD LINK COMPLEMENT 

The Whitehead link complement is known to be complete hyperbolic [14]. This 
example considers Dehn surgery on only one of the two cusps of the Whitehead 
link complement. 
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QUANTITIES FROM SNAP 

n = 4 

h = 1 and k = 2 
a — 

[0.9999343700073827649570992430+1.000170536257729817727630077*1, 
0.4999147436597508540443693049+0.4999671844066970777583211769*1, 

0.5000852675298210651958243937+0.5000328032070212542658981140*1, 
0.4999147436597508540443693049+0.4999671844066970777583211769*1] 

F = 

[20983, 0, -9872, 0, -9872, 11111, -1239, 20983, -2; 
0, 0, 0, 1, 1, -1, 0, 0, 0] 

G = 

[1. 1. 1. 1. 1, -2, 0, 0, -1; 0. -1. -1. -1. -1. 1. 1. 1. 1; 
-1, 1, 1, 1, 1, 0, -2, -2, -1;0, -1,-1,-1, -1, 1, 1, 1, 1] 

Pari-Gp CALCULATIONS 

\b\ = 6.290546043622649509854067366063508951285^-24 
L = 56237.01131396100111291495604741250466464 
\f'{a)-^\,up = 1.063909899076773471157618529051471308315 
|/'(a)"^|ien = 1.235415661324873497175222236812823735348 

, = 0.000007854853193291278165225494981053686965848 

, = 0.000005825343870778317976532920417278552662252 

W lien-'" 

3. LARGELINK COMPLEMENT 

This is the smaller of two extremely large link complements. See figure 8. It 
has 32 tetrahedra and 4 cusps. The other one has 57 tctrahcdra and 11 cusps. 
These two links are used by Leininger [10] to construct other knots and links by 
cut and paste methods, and then looking at their covers. For any even integer 
> 0, we eventually get from Largelink a two component link whose complement 
in contains an embedded totally geodesic surface of genus g. The importance 
of Largelink is that prior to this, such embedded surfaces could only be found in 
the complement of links with more than two components. 

QUANTITIES FROM SNAP 

n = 32 
h = k = A 
a — 

[5.431680776271168985E-77+1.043190149785894973378994944*I, 
0.4788708557877967957032308372+0.4995533597773714501527030266*1, 

-4.471822153042346518E- 77+0.9585980084313877504633692171*1, 
0.5211291442122032042967691627+0.4995533597773714501527030266*1, 
0.2929970420861826752219808548+1.473911044296957810855392169*1, 
-0.4509782171525463654321193064+1.200765444220459728291241593*1, 
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Figure 8: The Link Largelink 



1.000000000000000000000000000+0.9300613056344272435239940348*1, 

0.4638110047891777790136229363+0.4986886369525889712130902195*1, 

0.2371128008078259554449875702+0.6313317290357266968810549581*1, 

0.3060049499572359024927254903+0.4485055715132523515850447879*1, 

0.8585412796265143611585133046+1.027932770073775455116265474*1, 

0.5000000000000000000000000000+0.4792990042156938752316846085*1, 

0.4375240155821198504166057790+0.8813536566549109733830907053*1, 

1.042258288424406408593538325+0.9991067195547429003054060533*1, 

0.6696776343174901312972923995+0.7426519144895069642064793083*1, 

0.4518888703362351094400929330+0.9102903934144040876554144906*1, 

0.03927009472823897821842546946+1.571359648665194056162768058*1, 

-8.46343996E- 78+0.9585980084313877504633692171*1, 

1.000000000000000000000000000+1.043190149785894973378994944*1, 

0.02451089142372681728394034675+0.5982980953722294364585245396*1, 

0.9621628947892310086730291057+0.3083453854492406606721071067*1, 

0.7354295168083648566686302069+0.5515583107626382072967381105*1, 

0.6213864977872760582396031709+0.4161571993484503024065288682*1, 

0.6756917822944407062548472825+0.1978399260627268524593119332*1, 

0.5213536432299720005859050458+1.346701507985612627940863123*1, 

0.2659365860052524158474189000+0.5690611275237113909012011313*1, 

0.8916797222785394793793396465+0.5330292860478110834980601119*1, 

0.4489838724616496515202858898+0.4713823217067450930880172825*1, 

0.5364433482241135276673629307+0.6234802797569418514001720639*1, 

-0.4489884234808609710528328543+0.3884305318039174460267193001*1, 

1.000000000000000000000000000+0.9585980084313877504633692171*1, 

-8.49494342E- 77+1. 043190149785894973378994944*1] 



[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 

-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1; 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 
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0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, -1, -1, 1, -1, -1, 0, 1, 1, 1, 2, 0, 1, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 

1, 0, -2, 0, 0, -1, 0, 0, 0, 0, 1, 1, -1, 1] 

-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1; 0, 0, 0, 0, 0. 0. 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 
0, 0, 0, 0, 0, 0, -1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0; -1, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, -1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0; 0, -1, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1; 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1; 0, 1, -1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, -1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 1, 0, -1, 1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 2; 0, 0, 1, -1, 

0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, -1, 0, 
0, -1, 0, 0, 0, 0, 1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0; 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, -1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, -1; 0, 0, 0, 0, 1, 1, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, -2; 0, 0, 0, 0, 0, -1, -1, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 2: 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 1, 0, -1, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, -1; 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1; 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 

0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 1; 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, -1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, -1, 0, 0, 0, -2; 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 

1, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, ^1, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, -1, 0, 
0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 1, 0, 1, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 

1, 0, 0, -1, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 1, -1, 0, 0, 0; 0, 
0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, -1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 1, 1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1; 0, 0, 0, 0, 
0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, -1, 0, 0, 
0, 0, 1, 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 
0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0; 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, -1, -1, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, -1, 0, 1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, -1, 1, 0, 0, 0, 1, 0, 0, 2; 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 1, -1, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 
0, 0, 0, -1] 

Pari-Gp CALCULATIONS 

\b\ = 2.890741236697218507543429035402903716418£;- 27 
L = 38.46960927036768465200292167581178343887 
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sup = 8.212846275527759925085525656342053316915 



1 



i|ien = 10.32145710779244812406937753131330598443 

0.0001926925132239904423664849871566682428236 



1 



0.0001220029142841818172845137711227723107218 



4.3 Cusped Census 

We can apply the tests of Theorem 1.2 to every manifold in the SnapPea cusped census. 
The results are found in the following theorem. 

Theorem 4.1 Every manifold in the SnapPea cusped census has a complete hyperbolic 
structure. 

A program was written in Perl [1] that issues commands to Snap to send tetrahedron 
shapes and filling equations for each manifold in the cusped census to an output file. 
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Then a Pari-Gp program reads the file, getting the needed data per manifold, and 

applies the template using this input. The program then prints out the results. The 
first run of this process determined that all but four manifolds, 5 168, 6 297, 7 1431 
and 7 1927, have a complete hyperbolic structure. The program rejected these four 
because each one, upon triangulation by Snap, had one tetrahedron shape parameter 
with an imaginary component that was effectively zero. This was remedied by revising 
the original Perl program to process only these four manifolds, and including the "ran- 
domize" command to get a different, acceptable triangulation. The Pari-Gp program, 
also revised to process only these four manifolds, was then run using the second Perl 
output file. The result was a determination that they also have a complete hyperbolic 
structure. 

These programs can be adapted to give other information, such as the maximum 
value that normb, the norm of b, assumes over all the manifolds in the cusped census. 
Call this maxnormb. Similarly, for each manifold, we can ascertain the larger of the 
two values that normb is compared to, and then the minimum of these maximum 
comparison values over all the manifolds in this census. We do this because as long 
as normb of a manifold is less than the larger of the two comparison values for that 
manifold, the manifold will have a complete hyperbolic structure. Then if normb of a 
manifold in the census is less than the smallest of these maximum comparison values 
over the whole census, that manifold is guaranteed to have a complete hyperbolic 
structure. Gall this minimum of maximum comparison values minmaxvalue. It tells 
us the precision needed to evaluate a manifold in the census. We have 
maxnormb = 1.717844093022015223183888589087321425164875899778 E-26 
minmaxvalue = 0.00000147831677691814063380907736140260722549837777747014. 
Thus, the approximate solution given by SnapPea, which is given to 10 digits but 
is computed to an internal precision of at least 15 significant digits, is sufficient for 
use as our ai, . . . ,a„. It is interesting to see that the largest normb is considerably 
smaller than the smallest comparison value over the entire cusped census. The Perl 
programs and output files, as well as the Pari-Gp programs and log files, can be found 
at [12]. These Perl programs also include data with respect to a third test for a solution 
to the equations. However, the third test yields smaller comparison values than the 
Kantorovich tests, so it has no effect on minmaxvalue. 
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